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Abstract 



We consider Higgs boson production through gluon-gluon fusion in hadron collisions, 
when a veto is applied on the transverse momenta of the accompanying hard jets. We 
compute the QCD radiative corrections to this process at NLO and NNLO. The NLO 
calculation is complete. The NNLO calculation uses the recently evaluated NNLO soft 
and virtual QCD contributions to the inclusive cross section. We find that the jet veto 
reduces the impact of the NLO and NNLO contributions, the reduction being more 
sizeable at the LHC than at the Tevatron. 
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1 Introduction 



The mechanism of electroweak symmetry breaking is one of the main issues in high-energy particle 
physics. In the Standard Model (SM) and in its supersymmetric extensions the mechanism is 
accomplished by elementary scalar doublets. The Higgs boson(s) |l| are thus a fundamental 
ingredient of the theory, and their search is the main goal of current and future colliders. 

The direct search at LEP implies a lower limit of Mh > 114.1 GeV (at 95% CL) || on the mass 
Mh of the SM Higgs boson, and shows an excess of events, which may indicate the production 
of a Higgs boson with mass near 115 GeV 0,^,0. Global SM fits to electroweak precision 
measurements favour a light Higgs (Mh ^200 GeV) ||. 

After the end of the LEP era, the search for the Higgs boson will be carried out at hadron 
colliders. Depending on the luminosity delivered to the CDF and DO detectors during the forth- 
coming Run II, the Tevatron experiments can yield evidence for a Higgs boson with Mh ^ 180 GeV 
and may be able to discover (at the 5<r level) a Higgs boson with M H ^ 130 GeV ||. At the LHC, 
the SM Higgs boson can be discovered over the full mass range up to Mh ~ 1 TeV after a few 
years of running . 

The dominant mechanism for SM Higgs boson production at hadron colliders is gluon-gluon 
fusion through a heavy-quark (top-quark) loop |§. This mechanism is often called direct Higgs 
production, to distinguish it from associated production of Higgs boson and vector bosons, heavy 
quarks, jets, and so forth. 

At the LHC |§, gg fusion exceeds all the other production channels by a factor decreasing 
from 8 to 5 when Mh increases from 100 to 200 GeV. When Mh approaches 1 TeV, gg fusion still 
provides about 50% of the total production cross section. At low values of Mh the dominant decay 
mode is H — ► bb, but it is swamped by the QCD background. This decay mode is thus exploited 
in the low-mass range 100 GeV ^ Mh ^ 120 GeV only in the case of the associated production 
Hti. In the case of direct production, the most important decay channels for Higgs searches |7j 
are the rare decay H — > 77 in the low-mass range 100 GeV ^ Mh ^ 140 GeV, and the channel 
H — ► ZZ — > 41 in the mass range 130 GeV ^ Mh ^ 700 GeV. In the intermediate-mass region 
150 GeV &M H ^ 190 GeV, direct Higgs production followed by the decay H — > WW —> 
is also relevant. The strong angular correlations of the final-state leptons from W decay are an 
important ingredient for this discovery channel |T0|] . 

At the Tevatron, gg fusion remains the main production channel f6j, but it is experimentally 
less important than at the LHC because the decay rate H — > 77 is too low to be observed. When 
Mh ^5 135 GeV, the most promising discovery mechanism |J is thus the associated production 
(qq — » V* — > HV) of the Higgs boson with a vector boson V (V = W or Z), whose leptonic 
decay provides the trigger for the signal from H — > bb. Nevertheless, since the decays H — > 
W*W*,Z*Z* have increasing branching fractions for M H ^130 GeV, gg -> H -> W*W*,Z*Z* 
are the natural channels to consider when 140 GeV ^ Mh ^ 180 GeV. In particular, the decay 
mode W*W* — > l + l~vv is quite important 0, since it is cleaner than W*W* — > Ivjj, and the 
decay rate H — > W*W* is higher than H — > Z*Z* by about one order of magnitude. 

An important background for the direct Higgs signal H — > W*W* — > l + l~vv is tt production 
(tW production is also important at the LHC), where t — > Wb, thus leading to b jets with high 
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in the final state. If the b quarks are not identified, a veto cut on the transverse momenta of the 
jets accompanying the final-state leptons can be applied to enhance the signal/background ratio. 
Imposing a jet veto turns out to be essential, both at the Tevatron || [□]] and at the LHC @, to 
cut the hard b jets arising from this background process. 



The next-to-leading order (NLO) QCD corrections to gg fusion are large [12|, [13, Approx- 
imate evaluations JTH] of higher-order terms suggest that their effect can still be sizeable. The 



computation of the next-to-next-to-leading order (NNLO) corrections is thus important to better 
estimate the cross section for direct Higgs production. 



Recently two groups |16j [17| have performed a first step in this direction, by evaluating the 
soft and virtual contributions to the NNLO partonic cross section a(gg — > H + X) in the large- 
M top approximation. Our calculation |16| was performed by combining Harlander's result 
for the two-loop amplitude gg — > H with the soft factorization formulae for tree-level UlJl 
{gg — > Hgg, Hqq) and one- loop |21 
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Hg) amplitudes, and then using the technique of 
Ref. p3]j . The independent calculation of Ref. used a different method, and the analytical 
results fully agree. 



In Ref. |16| we also studied the quantitative impact of the soft and virtual NNLO contributions 
on direct Higgs boson production at the LHC. This was consistently done by using the recent 
MRST2000 set of parton distribution functions [53], which includes (approximated NNLO 



densities. In this paper we first perform an analogous study[] at the Tevatron Run II, and we show 
that the soft and virtual contributions are expected to give a fairly good estimate of the complete 
NNLO result. 

In the second part of the paper we study the effect of a jet veto on direct Higgs production 
beyond the leading order (LO). We present new NLO and NNLO calculations for the vetoed cross 
section both at the Tevatron and at the LHC. These calculations are performed by subtracting 
the LO (NLO) cross section for the production of Higgs plus jets from the corresponding inclusive 
NLO (NNLO) cross section. The LO and NLO subtracted cross section for the production of Higgs 
plus jets is evaluated in the large- M top limit, but without any soft and virtual approximations. At 
LO we derive analytical expressions, while at the NLO we use the numerical program of Ref. [J27[| . 

The paper is organized as follows. In Sect. |2| we recall the theoretical framework and the 
analytical results for the soft and virtual NNLO corrections. In Sect. [3] we present the results for 
inclusive Higgs production at the Tevatron Run II. In Sect. |] we present the results obtained by 
applying a jet veto on the inclusive cross section, both at the Tevatron Run II and at the LHC. 
Our conclusions are summarized in Sect. |5|. 



tPart of these results was anticipated in Ref. 
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2 Inclusive QCD cross section at NNLO 



We consider the collision of two hadrons hi and h 2 with centre-of-mass energy sfs. The inclusive 
cross section for the production of the SM Higgs boson can be written as 



a ,b Jo 



fb/h 2 (x 2 , Hp) I dz 61 z 



XiX 2 



a z G ab (z; a s (fi 2 R ), M 2 H / fi 2 R ; M 2 H /fi 



F) > 



(1) 



where th = M^/s, and fip and ixr are factorization and renormalization scales, respectively. The 
parton densities of the colliding hadrons are denoted by fa/h^iHp) an d the subscript a labels the 
type of massless partons (a = g, qf, q~f, with Nf different flavours of light quarks). We use parton 
densities as defined in the MS factorization scheme. 

From Eq. ([!]) the cross section a a b for the partonic subprocess ab — * H+X at the centre-of-mass 
energy s = Xix 2 s = Mjj/z is 



a ab (s,M 2 H ) = — a M 2 H G ab (z) = a z G ab (z) 

s 



(2) 



where the term 1/s corresponds to the flux factor and leads to an overall z factor. The Born-level 
cross section o"o and the hard coefficient function G ab arise from the phase-space integral of the 
matrix elements squared. 



Equation (p]) can also be recast in the form 

dr 



/ dr 
— £ab/hih 2 ( T , /4) G ab (T H /r) 



where the function 



1 dx 

Cab/hihair, fip) = I —f a/hl (x,fl 2 F )f b/h2 (T/x,fl 2 F ) 



(3) 



(4) 



is the parton luminosity that weights the event initiated by partons a and b. 



The incoming partons a, b couple to the Higgs boson through heavy-quark loops and, therefore, 
<7 and G ab also depend on the masses Mq of the heavy quarks. The Born-level contribution cr 
is IB 
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where Gp = 1.16639 x 10 5 GeV 2 is the Fermi constant, and the amplitude Aq is given by 



Aq(x) = -x 



l + (l-x)f(x) 

. n 1 



arcsm 



in , — tir 

1 - V^x 



x>l 



x < 1 



(6) 
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In the following we limit ourselves to considering the case of a single heavy quark, the top quark, 
and Nf = 5 light-quark flavours. We always use M top (M top = 176 GeV) to denote the on-shell 
pole mass of the top quark. 

The coefficient function G ab in Eq. (|l]) is computable in QCD perturbation theory according 
to the expansion 



G ab (z;a s (^ 2 R ),M 2 H /fi 2 R ;M 2 H /^ 



a, 



g|04) G w 



7T 



ab 



i(/4)E 
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_ M 2 H _ M 2 H 
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7T 
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«s(/4) ^(2) 



7T^ 



/4 /4 



C»(a 5 



where the (scale-independent) LO contribution is 

G ( ' b \z)=5 ag 5 bg 5(l-z). 



(7) 



(8) 



Throughout the paper we work in the framework of the large-M top expansion, where one 
can exploit the effective-lagrangian approach |28|, [ 



15| to embody the heavy-quark loop in an 



effective vertex. However, unless otherwise stated, we include in a the full dependence on M top . 
This approximation |14], IE turns out to be very good when M H < 2M top , and it is still accurate]] 
to better than 10% when M H &1 TeV. 



The NLO coefficients G^ b \z) are known [0, [14|]. Their explicit expressions (see e.g. Ref. JT 
contain three kinds of contributions: 



Virtual and soft contributions, which are respectively proportional to 5(1 — z) and to the 
distributions T>Az), where 



V, 



ln l (l 



These are the most singular terms when z 



l-z 
-> I. 



(9) 



Purely collinear logarithmic contributions of the type ln l (l — z). These contributions give 
the next-to-dominant singular terms when z — > 1. 

Hard contributions, which are finite in the limit z — > 1. 



The soft- virtual (SV) approximation is defined ]T6[ by keeping only the terms proportional to 
5(1 — z) and T>i(z) in the coefficient G ab . In this approximation only the gg channel contributes 
and we have 



Gab SV ( Z 'i M h/^r'i M h/Vf) = $agh 



x Al x 33-2iV f , u 2 

5(1 -z) (y + 6C(2) + ^^m| 

M 2 



+6V In— £ + 12 V 1 



(10) 



*The accuracy of this approximation when Mh & 2M top may not be accidental. In fact, as discussed in the 
following, the main part of the QCD corrections to direct Higgs production is due to parton radiation at relatively 
low transverse momenta. Such radiation is weakly sensitive to the mass of the heavy quark in the loop. 
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(2) 

The NNLO coefficients G ab are not yet completely known, but their SV approximation has 



been computed in Refs. |16|, [L7j • It reads: 
Gf^(zM/^ R} M] I /^ F )=6(l-z) 



+ V 
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(11) 



To give a rough idea of the numerical hierarchy of the LO, NLO and NNLO contributions, we 
can set fip = = M H = 115 GeV and as(M H ) = 0.112. The first three terms in the expansion 
(0) thus give 



(4^6(1 -z) ( 1 + 0.548 + 0.105 

+v ( o + o + 0.283 ; 

+Vi ( + 0.428 - 0.040 

+v 2 ( o + o - 0.029 ; 

+V 3 ( + + 0.092 ' 



Since 1 > z = th /t > t#, the SV terms are certainly the dominant contributions to the cross 
section in the kinematic region near the hadronic threshold (t# = M\j s ~ 1). At fixed s, this 
means that these terms certainly dominate in the case of heavy Higgs bosons. However, the SV 
terms can give the dominant effect even long before the threshold region in the hadronic cross 
section is actually approached P0f . This is because, in the evaluation of the hadronic cross section 
in Eq. (|3]), the partonic cross section a a b(s, Mjj) has to be weighted with the parton luminosities, 
which are strongly suppressed at large r = s/s = X\X 2 - In other words, owing to the strong 
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suppression of the parton densities f a (x,fXp) at large x, the partonic centre-of-mass energy a/I is 
typically substantially smaller than yfs ((s) = (xix 2 s) = (r)s), and the variable z = M\js in 
G a b(z) can be close to unity also when y/s is not very close to M H . 



10 3 r S .M H =150 GeV MRST2000 LO - 




200 400 

Vs [GeV] 



Figure 1: Parton (gg and qg+gq) luminosities jC(t, fi F ) with = Mh = 150 GeV at the Tevatron 
Run II as a function of the partonic centre-of-mass energy \/§ (s = rs). We use the LO set of 
parton distributions in Ref. |||/. 



At fixed Mh, the quantitative reliability of the large- ,2 approximation depends on the value 
of \/i and on the actual value of the parton luminosities C a b/h 1 h 2 ( T , A*f) i n Eq. @. In Fig. [I] the 
gg and (qg + gq) luminosities (with fip — Mh = 150 GeV) at the Tevatron Run II are plotted 
as a function of the partonic centre-of-mass energy \f§ (s = rs). We see that the luminosities 
decrease by almost two orders of magnitude when a/I increases from Mh = 150 GeV (hadronic 
threshold) to 300 GeV. A similar effect occurs at the LHC (see Fig. 1 in Ref. [2B]). We thus expect 



that the large-z expansion of the coefficient function G a b(z) reliably approximates the size of the 
higher-order QCD corrections to Higgs boson production at the Tevatron and the LHC. At NLO, 



this was explicitly checked in Ref. fL5f , by comparison with the complete expression of G 



ab 



The authors of Ref. |I5] also pointed out that at NLO the numerical effect of the logarithmic 
term ln(l — z) of collinear origin is not small. Following this observation, in Ref. |16| we intro- 
duced the soft-virtual-collinear (SVC) approximation of the coefficient function G a b{z). The SVC 
approximation is defined by including the leading ln fc (l — z) contribution from the collinear region 
in the gg channel: 

G^ c (z; Ml, 74; Mil '^ F ) = G^(z; M 2 H /^ R ; M 2 H /» 2 F ) - 12 ln(l - z) , (12) 



Gf^ S VC (z; M|//4 MH A) = G(f v (z; M 2 H /fa Ml/u 2 ,) ~ 72 ln 3 (l - z) . (13) 

Since the term ln fc (l — z) added to the SV expressions is that with the highest power k at each 
perturbative order (k = 1 and k = 3 at LO and NLO, respectively), the SVC approximation 

(2) 

consistently includes the next-to-dominant contribution to G gg in the limit z — > 1. The comparison 
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between the SV and the SVC approximations can thus be used |T6| to gauge the quantitative 
accuracy of approximating G a b(z) by its large- z limit. In the next section we study the impact of 
the SV and SVC approximations at the Tevatron Run II. 

Two different large-z approximations, named 'soft' and 'soft+sl' were also considered in the 
numerical study of Ref. |17|. The 'soft' approximation of Ref. fl?]] regards the whole partonic cross 



section a a b in Eq. (0), while our SV approximation refers only to the hard coefficient function 
G a b{z). In other words, we perform the soft approximation on the phase-space integral of the 
matrix elements squared, while the kinematical flux factor 1/s in a a b is kept fixed and not expanded 



around s = s. This means that we expand G a b{z) around z = 1, whereas the authors of Ref. |L7 
consider the expansion of G a b(z) = zG a b(z). Owing to the identity 

*X\ = Z\-ln*(l-z) , (14) 

the two expansions differ by subdominant ln l (l — z) contributions of kinematical origin (see also 
Sect. 3 in Ref. jl6|). Moreover, the 'soft+sl' approximation of Ggg |17j includes also additional 



subleading terms, proportional to ln 2 (l — z) and ln(l — z), whose coefficients (unlike those in 
Eqs. ( |ITD and are not (exactly) predictable at present | I5| . In summary, we note that there 



is no one-to-one correspondence between 'soft' (or 'soft+sl') in Ref. fl7| and SV (or SVC) in our 
definition, the difference being due to logarithmic contributions that, formally, are consistently 
subdominant when z — > 1. 

Although, from a formal viewpoint, it is legitimate to perform the large- z expansion of either 
G a b( z ) or G a b(z), the two expansions can quantitatively differ when applied to the evaluation of 
the hadronic cross section in Eq. (|l|). In fact, when using G a b(z), the analogue of Eq. (|j) is 

/ 1 dr ~ ~ 
£ab/h 1 h 2 (l',IJ. 2 F )G ab (TH/T) . (15) 

Comparing Eqs. (|3|) and ([[§), we see that the coefficient function G a b{z = t h /t) is convoluted 
with the momentum-fraction luminosity C(t, up) = tC(t,/j, f ). Owing to the relative rescaling 
factor t = s/s, £(r, fi F ) is much less steep than the luminosity £(r, ji 2 F ) (see Fig. []]) that enters 
Eq. (H]). Therefore, in the numerical evaluation of the Higgs boson cross section, we expect and 



anticipate (see the comment at the end of Sect. |) that the large-z expansion of G a b(z) |17 
converges more slowly than that of G a b(z). The slowing down is ultimately caused by a too 
extreme kinematics approximation of the partonic cross section in Eq. (0): the flux factor 1/s has 
been replaced by 1/Mjj. 



3 Inclusive production at the Tevatron Run II 

In this section we study the higher-order QCD corrections to the inclusive production of the SM 
Higgs boson at the Tevatron Run II, i.e. proton-antiproton collisions at ^/s = 2 TeV. We recall 
that we include the exact dependence on M top in the Born- level cross section o"o (see Eq. ([5])), 
while the coefficient function G a b(z) is evaluated in the large-M top limit. At NFO JTJ], |T5| this is 
a very good numerical approximation when Mh < 2M top . 
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Unless otherwise stated, cross sections are computed using the MRST2000 [^4| sets of parton 
distributions with densities and coupling constant evaluated at each corresponding order, i.e. using 
LO distributions and 1-loop «s f° r the LO cross section, and so forth. The corresponding values 
of J^qcd [a a (M z )) are 0.132 (0.1253), 0.22 (0.1175) and 0.187 (0.1161) GeV, at 1-loop, 2-loop 
and 3-loop order, respectively. In the NNLO case we use the 'central' set of MRST2000, obtained 
from a global fit of data (deep inelastic scattering, Drell-Yan production and jet Et distribution) 
by using the approximate NNLO evolution kernels presented in Ref. p5[ . The result we refer to 
as NNLO-SV(SVC) corresponds to the sum of the LO and exact NLO (including the qg and qq 
channels) contributions plus the SV(SVC) corrections at NNLO, given in Eq. (|TT|) (Eq. (|T3|)). 
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Figure 2: Relative difference of LO and NLO results computed with GRV98 and CTEQ5 distribu- 
tions with respect to MRST2000 distribution. 



In Fig. we show the dependence of the LO and NLO cross sections on the choice of parton 
distributions. The results obtained by using both the GRV98 |31[| and CTEQ5 |3^] sets differ 



substantially from the MRST2000 results. In the case of the GRV98 sets, this difference is not 
unexpected, since the value of as(Mz) and the gluon distribution are quite different from those of 
MRST2000 and CTEQ5. We see that, as M H increases from 100 to 200 GeV, the LO result ob- 
tained by using the CTEQ5 set is between 10 and 30% lower than the one obtained by MRST2000. 
At NLO the difference is smaller, and it is compatible with the ±10% uncertainty recommended by 
the CTEQ collaboration [[^J on gg and qg luminosities in the x region that controls Higgs boson 
production at the Tevatron. Although the LO differences are physically less meaningful than the 
NLO differences, the size of the former has to be kept in mind when quoting QCD predictions 
(based on either analytic calculations or Monte Carlo event generators) that directly or indirectly 
(e.g. in the case of ^-factors) use LO parton densities. 

Figure ^ shows the scale dependence of the cross section for the production of a Higgs boson 
with Mh = 150 GeV. The scale dependence is analysed by varying the factorization and renor- 
malization scales by a factor of 4 up and down from the central value Mh- The plot on the left 
corresponds to the simultaneous variation of both scales, Hf = = X^h, whereas the plots in 
the centre and on the right respectively correspond to the results of the independent variation of 
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Figure 3: Scale dependence of the Higgs production cross section for M H = 150 GeV at LO, NLO, 
NNLO-SV and NNLO-SVC. 

the factorization and renormalization scales, keeping the other fixed at the central value. 

As expected from the QCD running of as, the cross sections typically decrease when fj, R 
increases around the characteristic hard scale Mh- A similar behaviour^] is observed when /ip 
varies, since the cross sections are mainly sensitive to partons with momentum fraction x ~ 0.05- 
0.1, and in this x range scaling violation of the parton densities is (slightly) negative. The largest 
variations in the cross section calculation are thus obtained by simultaneously varying and 
fip- The scale dependence is mostly driven by the renormalization scale, because the lowest-order 
contribution to the process is proportional to «g, a (relatively) high power of as- 

In summary, Fig. |3] shows that the scale dependence is reduced when higher-order corrections 
are included. Varying the scales in the range 0.5 < \Ri Xf < 2, the reduction is from about ±20% 
at full NLO to about ±10% and ±15% at NNLO-SV and NNLO-SVC, respectively. The increase 
in the scale dependence when going from NNLO-SV to NNLO-SVC is due to the fact that the 
contribution of the dominant collinear terms included in the SVC approximation (see Eqs. (0) 
and (|T3"|)) is not small and scale-independent, so it cannot be compensated by scale variations. 

In Fig. |] we study the .fT-factors, defined as the ratio of the cross section evaluated at each 
corresponding order over the LO result. Since the LO result sizeably depends on the choice 
of parton distributions, the i^-factor should be interpreted with care. For instance, for Mh = 
180 GeV, the NLO if-factor computed by using MRST distributions is about 20% smaller than 
the one obtained by using CTEQ distributions. 

In Fig. f| the bands account for the 'theoretical uncertainty' due to the scale dependence, 
quantified by using the minimum and maximum values of the cross sections when the scales /ar 
and [if are varied (simultaneously and independently, as in Fig. [D in the range 0.5 < XRi Xf < 2. 
The LO result that normalizes the .^-factors is computed at the default scale Mh in all cases. 

§At the LHC, the \if dependence is opposite, as shown in Fig. 1 of Ref. p6| . 



9 




Figure 4: K-factors for Higgs production for the full NLO result and the NLO-SV, NLO-SVC, 
NNLO-SV and NNLO-SVC approximations. 

The plot on the left-hand side of Fig. |] shows the scale uncertainty at LO and compares the 
full NLO result with the NLO-SV and NLO-SVC approximations. When Mh is in the range 
140-180 GeV, the NLO-SV approximation tends to underestimate the full NLO result by about 
7-8%, whereas the NLO-SVC approximation overestimates it by about 4-5%, showing the effect 
of the term ln(l — z) added in the SVC approximation. We find that the SVC result provides 
an excellent approximation of the full gg contribution at NLO, the difference being less than 2 
per mille. The contribution of the qg channel tends to lower the SVC result by ^5%. The effect 
of the qg channel is mainly due to the logarithmically-enhanced behaviour G^qg{z) ~ | ln(l — z) 
of the corresponding coefficient function at large z. The improved reliability of the SV and SVC 



approximations with respect to the LHC |16| is not unexpected since at the Tevatron we are closer 
to threshold. 

The right-hand side of Fig. ^ shows the SV and SVC results at NNLO. Again, the SVC band 
sits higher than the SV one and, as shown in the inset plot, the ratio of the corresponding cross 
sections increases with respect to NLO. The contribution from non-leading terms ln fc (l — z), with 
k < 3 (which are not under control within the SVC approximation), is not included. We have 
tried to add a term ln 2 (l — z) with a coefficient as large as the one of the ln 3 (l — z) contribution, 
and we have found a small ((9(5%)) modification. Therefore, we expect the effect of these non- 



leading logarithmic terms to be numerically less important. As pointed out in Ref. |16[, at NNLO 
there is a (still unknown) leading collinear contribution proportional to ln 3 (l — z) also in the qg 
channel. Owing to the size of the contribution of the qg channel at NLO, the quantitative effect 
of this NNLO term can be of the same order as the corresponding one included through the SVC 
approximation in the gg channel. 

We recall that the NNLO results on the right-hand side of Fig. ^ are obtained by using the 
(approximated) NNLO parton distributions of the MRST2000 set. Using the NLO parton distri- 
butions, the i^-factors would be smaller by about 5 to 8% as Mh increases from 120 to 180 GeV. 
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Mu(GeV) 


LO 


NLO 


NLO-SV 


NLO-SVC 


NNLO-SV 


NNLO-SVC 

lili J_J V. J kJ V v_. 


140 


2282 


0.4715 


0.4338 


0.4922 


0.6163 


0.7222 


150 


0.1856 


0.3794 


0.3507 


0.3969 


0.5002 


0.5841 


160 


0.1523 


0.3080 


0.2859 


0.3229 


0.4092 


0.4763 


170 


0.1260 


0.2519 


0.2349 


0.2646 


0.3369 


0.3912 


180 


0.1050 


0.2075 


0.1943 


0.2184 


0.2795 


0.3235 



Table 1: Cross sections in pb as a function of Mh- The calculation is performed by setting 
Vr = = Ma- 



in summary, considering the results obtained at NLO, we expect the full NNLO i^-factor to be 
between the SV and SVC bands. For reference, we give in Table 1 the central values of the cross 
section in the range M H = 140-180 GeV. In particular, for a light Higgs boson (M H ^200 GeV), 
this corresponds to an increase of about 50% with respect to the full NLO result, i.e. a factor 
of ~ 3 with respect to the LO result. Taking into account that the NLO result increases the LO 
cross section by a factor of about 2, our result shows that the convergence of the perturbative 



series is poorer at the Tevatron than at the LHC [16(]. This also implies that QCD contributions 



beyond NNLO can still be significant at the Tevatron. 

Our conclusion on the possible relevance of higher-order contributions is not in contradiction 
with the improved scale dependence of the NNLO calculation. As is well known, the customary 
procedure of varying the scales to estimate the theoretical uncertainty due to missing higher-order 
terms can only give a lower limit on the 'true' uncertainty. This is well demonstrated by the plot 
on the left-hand side of Fig. [|, which shows no overlap between the LO and NLO bands. Since 
the NLO and NNLO bands still tend to show no (or a marginal) overlap, their size cannot yet be 
regarded as a reliable estimate of the theoretical uncertainty. 



A comment on the numerical results presented in Ref. [17] is in order. The authors of Ref. [17 



use different parton distributions and, as explained in Sect. |2], there is no one-to-one correspondence 
between their 'soft' (or 'soft+sl') approximation and our SV (or SVC) approximation. Therefore, 
those results cannot directly be compared with ours (see also Ref. p4|). Certainly, they show a 



slower numerical convergence of the large- z expansion used in Ref. f[? |. As discussed at the end 



of Sect. 0, this is not unexpected because of significant subdominant effects of kinematical origin. 



4 Vetoed cross section 



Direct Higgs production followed by the decay H — > WW — > l + l~vu is an important channel to 
discover the SM Higgs boson in the intermediate-mass range 140 GeV ^ Mh & 190 GeV both at 
the Tevatron and at the LHC. Nevertheless, as recalled in Sect. [I], several experimental cuts have 
to be applied to discriminate the signal over the background. An important selection to enhance 
the statistical significance is a veto on the high transverse-momentum jets that accompany the 
production of the Higgs boson. Events with high transverse-momentum jets are excluded from 
the analysis g [7J, [TJ. 

In this section we study the effect of a jet veto on inclusive Higgs production. The events that 
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pass the veto selection are those with p i ^ t < p™ to , where pi 1 ? is the transverse momentum of any 
final-state jets. Jets are defined by a cone algorithm. In the perturbative calculation, jets are 
represented by a parton or a set of partons. The Higgs boson can be accompanied by one final- 
state parton at NLO, and by one or two final-state partons at NNLO. In the case of a single final- 
state parton with transverse momentum p 1T , the vetoed cross section is computed by imposing 
I Pit | < Pr to (i-e. we veto events with | pi t| > p T eto ). When there are two final-state partons with 
transverse momenta pi t and P2T, we consider their angular distance R\ 2 = — f^) 2 + {4>\ ~ 'P2) 2 
in the pseudorapidity-azimuth plane: if i? i2 > R, we impose the constraints |pir|, |P2t| < Pr to 
(i.e. we veto events with |pir| > Vr to or |P2t| > Vr to )'i if R12 < R, we combine the two partons 
in a single jet and we impose |pit + P2t| < Vr to (i- e - we veto events with |pit + P2t| > p™ to ). In 
all the numerical results presented in this section, the cone size R of the jets is fixed at the value 
R = 0.4. 



As in the case of the inclusive cross section, we evaluate the jet-vetoed cross section by using 
the large-M top limit. Studies on Higgs+jets production show that the QCD matrix elements in 
the large- M top limit [^, approximate very well those with the exact M top dependence |37|, |38fl , 



provided the jet transverse momenta and, hence, p^ to remain small with respect to M top . 

The vetoed cross section <7 veto (s, M#; p™ to , R) can be written as 

a veto ( S , M 2 H - p™ to , R) = a(s, M 2 H ) - Aa(s, M 2 H ; ^ eto , R) , (16) 

where a(s,Mfj) is the inclusive hadronic cross section in Eq. (P, and Act is the 'loss' in cross 
section due to the jet-veto procedure. The vetoed cross section is computable by a factorization 
formula analogous to Eq. ([!]), apart from the replacement of the coefficient function G a b{z) by 
the vetoed coefficient function G^f°(z; itt,R)- By analogy with Eq. ([TB]), we can write the vetoed 
coefficient function as 



GT(z; ttt, R) = G ab (z) - AG ab (z; 7T T , R) , (17) 

where, to simplify the notation (see Eq. (]20D), the dependence on p™ to is parametrized by the 
dimensionless variable ttt'- 

o„vcto rz 

^•*"°/M„) - fcfe ■ (18) 



At LO, the vetoed cross section is equal to the inclusive one, so the subtracted coefficient 
function AG a b vanishes. At higher perturbative orders, AG a b is computable according to the 
power-series expansion 

AG ai (z- ttt, R; a s (»l), M' H /ii 2 H - M 2 „/^ F ) = AG™ (z; ir T ) 9(1 - wr(z, p?°/M H )) 

+ d^ A G$L« T , R M.,MPl + 0(4). (19) 

The NLO contribution AG^ is scale- independent and R- independent, but it depends on p™ to . 
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We have explicitly evaluated it in the large- M top limit, and we find: 



*GW(z;« T ) = P„{z) In 



7T 



T 



11 a-z) 



1 - ji - 4 2 



7T; 





z 



TTrp 

22 



1 -7r| , 



1 - 7r| , 



(20) 



32 



27 



1 - 



7T. 



T 



l-7r| , 



where 



Al + (l-z) 



PM = 6 



1 — z z .„ . 
+ + z(l - z) 



1 - z 



(21) 



The evaluation of the NNLO contribution AG 



(2) 



cannot be easily performed in analytic form, 

£j numerically 

by using the program of Ref. |57j|. The program, which implements the matrix elements (in the 
large- M top limit) of Refs. [^, by using the subtraction-method procedure of Ref. f39fl , computes 



ab 

since the calculation depends on the details of the jet algorithm. We compute AG ab 



the QCD corrections to Higgs+jet(s) production up to order a|. For the purpose of evaluating 
the contribution Act to Eq. (0), we use the program to compute the Higgs+jet(s) cross section 
when the transverse momentum of the highest-pr jet is larger than p™ to . 

In the following we present both NLO and NNLO numerical results for the vetoed cross section 



<7 veto in Eq. (|16|). The results are obtained by using the parton distributions of the MRST2000 
set, as explained in Sect. [| The NLO calculation is exact: apart from using the large- M top limit, 
we do not perform any further approximations. At the NNLO, the contribution Act to Eq. (|16|) 
is again evaluated exactly, while to evaluate the contribution of the inclusive cross section we rely 
on our approximate estimate in Sect. |3] (see the corresponding LHC results in Sect. 4 of Ref. JT5|]) 
and, in particular, we use the NNLO-SVC result. Therefore, once the full NNLO result for the 
inclusive cross section is available, it can straightforwardly be used to 'correct' our NNLO estimate 
for the vetoed cross section. As stated above, in our numerical calculations we fixed the cone size 
of the jets to the value R = 0.4. The independence of the perturbative calculation first appears 
at the NNLO. In particular, the NNLO vetoed cross section decreases by increasing R. 



Note that the numerical program of Ref. |27) is a Monte Carlo code that evaluates Higgs+jet(s) 
production at the fully exclusive level. Therefore, it can be used through the subtraction procedure 
of Eq. flTEp to compute vetoed cross sections also when additional kinematical cuts (e.g. cuts on the 
rapidities of the jets, or asymmetric cuts on the jet transverse momenta) or different jet definitions 
are considered. 

We first present the vetoed cross section^] at the Tevatron Run II. In Fig. [5] we show the 
dependence of the NLO results on the Higgs mass for different values of p™ to (15, 20, 30 and 

^The numerical program of Ref. |27j implements the large- Mt op limit strictly, i.e. also the Born-level contribution 
(To in Eq. (^) is evaluated in the limit Mh /M top — * 0. For simplicity, the numerical results of this section implement 
the same approximation. The approximation used in Sect. || can be recovered by simply rescaling the absolute 
value of the cross sections by the overall factor <j /<j (MH/M tO p = 0). Of course, such a rescaling has no effects on 
the ratios Aa/a and on the vetoed K-i actors. 
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NLO MRST2000 ' : ha/a 
Mf=Mr=M h 4 " 



a[pb 




130 140 150 160 170 180 1£ 
M H (GeV) 



1.5 



0.5 



K— factor 



p T <50 GeV- 
p x <30 GeV- 

p T <20 GeV- 
p T <15 GeV- 



130 140 150 160 170 180 190 
M H (GeV) 



Figure 5: Vetoed cross section and K-factors: NLO results at the Tevatron Run II. 

50 GeV). The vetoed cross sections cr veto (s, Mfj]p^ to , R) and the inclusive cross section a(s,M]j) 
are given in the plot on the left-hand side. The inset plot gives an idea of the 'loss' in cross section 
once the veto is applied, by showing the ratio between the cross section difference Act in Eq. (|16|) 
and the inclusive cross section at the same perturbative order. As can be observed, for large 
values of the cut, say p^ to = 50 GeV, less than 10% of the inclusive cross section is vetoed. The 
veto effect increases by decreasing p™ to , but it is still smaller than 30% when p™ to = 15 GeV. On 
the right-hand side of Fig. |5|, we show the corresponding X-factors, i.e. the vetoed cross sections 
normalized to the LO result, which is independent of the value of the cut. 




130 140 



150 160 170 

M H (GeV) 



180 190 



K— factor 



Inclusive 



p T <50 GeV 

p x <30 GeV 

p T <20 GeV 

p T <15 GeV 



130 140 150 160 170 180 190 
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Figure 6: Vetoed cross section and K-factors: NNLO results at the Tevatron Run II. 



Figure |5] shows the analogous results at NNLO. Note that, although Aa/a is slightly larger 
than at NLO, the vetoed cross sections and, thus, the iv'-factors, are still larger than at NLO. This 
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is mostly due to the large increase of the inclusive cross section at NNLO, as shown in Sect. |[ 




Figure 7: Scale dependence of the inclusive and vetoed cross sections at LO, NLO and NNLO-SVC 
at the Tevatron. 

All the results plotted in Figs. [| and |6| have been obtained by fixing the renormalization and 
factorization scales at the default value = Hf = M H . In Fig. [7j, we study the scale dependence 
of the cross sections (inclusive and with p™ t0 = 15 and 30 GeV) at LO, NLO and NNLO. The 
renormalization and factorization scales are varied simultaneously by a factor of 4 up and down 
with respect to Mjj. No significant differences are found when the scales are varied separately. 
As can be observed, there is an improvement in the stability of the result when higher order 
corrections are included. This is particularly noticeable when going from LO to NLO, while the 
comparison between NNLO and NLO looks rather similar to the one observed in the inclusive 
case. 

The scale-dependence effects on the perturbative .fT-factors can be appreciated also from Fig. |8], 
where the LO, NLO and NNLO-SVC bands (computed as in Fig. ^) are shown in the case in which 
a veto of p^ to = 15 GeV is applied. Comparing Fig. |8] with the inclusive case in Fig. ^, we see that 
the effect of the veto is to partially reduce the relative difference between the NLO and NNLO 
results; the increase of the corresponding i^-factors can be estimated to about 25%. 

The results on the jet veto presented so far can be qualitatively explained by a simple physical 
picture. As shown in Sect. |||, the effect of the higher-order contributions to inclusive Higgs boson 
production at the Tevatron is large. The dominant part of this effect is due to soft and collinear 
radiation, whereas the accompanying hard radiation has little effect. The characteristic scale of 
the highest transverse momentum p™ ax of the accompanying jets is indeed p™ ax ~ (1 — z)Mh 
(see e.g. the upper bound on p^ to from the theta function in Eq. (|l9l)), where the average value 
(1 — z) — (1 — Mjj/s) of the distance from the partonic threshold is small. As a consequence the 
jet veto procedure is weakly effective unless the value of p™ to is pretty small (i.e. substantially 
smaller than p™ ax ). Decreasing p™ to , the enhancement of the inclusive cross section due to soft 
radiation at higher orders is reduced, and the jet veto procedure tends to improve the convergence 
of the perturbative series (see e.g. Figs. [| and Note also that the characteristic scale 
is a slightly increasing function of Mjj, the linear increase with Mh being partially compensated 
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Fi gure 8: K -factors for Higgs production at the Tevatron in the case of a veto of p™ to 
LO, NLO and NNLO-SVC approximation. 



15 GeV: 



by the decrease of (1 — z). Therefore, at fixed p™ to the vetoed .fT-factor decreases more than the 
inclusive i^-factor when Mh increases (see Fig. 

On the basis of this physical picture, we can easily anticipate the qualitative effect of the jet 
veto at the LHC The overall features of the QCD corrections to inclusive Higgs boson production 



at the LHC []I6[ are the same as at the Tevatron. The main quantitative differences arise from 
the fact that Higgs production at the LHC is less close to threshold than at the Tevatron and, 
therefore, the accompanying jets are less soft ((1 — z) is larger) at the LHC than at the Tevatron. 
At fixed p™ to the effect of the jet veto is thus stronger at the LHC. 
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Figure 9: Vetoed cross sections and K-factors at NLO at the LHC. 
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Figure 10: Vetoed cross sections and K-factors at NNLO at the LHC. 

The results for the vetoed cross sections at the LHC are presented in Figs. |^ and [It] for 
p™to _ 20 ; 30, 50 and 70 GeV. At fixed value of the cut, the impact of the jet veto, both in 
the 'loss' of cross section and in the reduction of the i^-factors, is larger at the LHC than at the 



Tevatron Run II. This effect can also be appreciated by comparing Fig. |TT] and Fig. [8]. At the 
LHC, the value of p^ to = 30 GeV is already sufficient to reduce the difference between the NNLO 
and NLO results to less than 10%. 
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Figure 11: The same as in Fig. ]Q, but at the LHC and with p^r = 30 GeV. 



Note that when p™ to is much smaller than the characteristic scale p™ ax ~ (1 — z)M H , the 
coefficients of the perturbative expansion of the vetoed cross section contain large logarithmically- 
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enhanced contributions. For instance, from Eq. (]20Q we have 



AG g 1 g \z;TC T )~2P gg (z) In 



(1 - z)M H 



^vcto 
FT 



(22) 



The presence of these contributions can spoil the quantitative convergence of the fixed-order 
expansion in as- Since (1 — z)Mh is larger at the LHC than at the Tevatron, the value of p™ to 
at which these effects become visible is larger at the LHC. The perturbative if-factors shown in 
Fig. [12| suggest that at the LHC the vetoed cross section is sensitive to these large logarithmic 
terms already when p™ to = 15 GeV. Indeed, the scale-dependence band is larger at NNLO than 
at NLO. At such small values of p^ to , perturbative contributions beyond the NNLO can still be 
significant. 
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Figure 12: The same as in Fig. ]Q, but at the LHC and with = 15 GeV. 



5 Conclusions 



In this paper we have presented our study of QCD radiative corrections to direct Higgs boson 
production at hadron colliders and the impact of jet veto. 



Using the theoretical results of Refs. [16|, |T^], in Sect. |3] we have first shown that the NNLO-SV 
and NNLO-SVC approximations are expected to be a good estimate of the full NNLO contributions 
to the inclusive cross section at the Tevatron Run II. Similar conclusions about the inclusive 



production at the LHC were obtained in Ref. [16]. The reliability of these NNLO approximations 



follows from the fact that a sizeable part of the higher-order QCD corrections to the partonic cross 
section is due to the emission of soft radiation, whereas the effect of hard radiation is suppressed by 
the steeply falling behaviour of the parton distributions at large x. As is customary practice, the 
perturbative expansion has been performed in the MS factorization scheme. Since the MS-scheme 
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parton densities embody (by definition) too much soft-gluon suppression (see e.g. Ref. [(Wj), the 
perturbative corrections to the partonic cross section have to compensate for it; they thus enhance 
the hadronic cross section. At the Tevatron, we have shown that the NNLO effect is large and 
increases the cross section by about 50% with respect to the NLO result. This suggest that 
contributions beyond the NNLO can still be significant. Since Higgs boson production is less close 
to threshold at the LHC than at the Tevatron, the accompanying QCD radiation is less soft at 



the LHC than at the Tevatron. The estimated size of the NNLO effects at the LHC |16] is thus 
consistently smaller than at the Tevatron. 

We have then addressed the impact of a jet veto on the inclusive cross section, both at the 
Tevatron and at the LHC. We have presented results at NLO and NNLO. At NNLO, we have used 
our estimate for the inclusive cross section and we have computed (subtracted) the effect of the jet 



veto by using the numerical program of Ref. [27]. The jet veto reduces not only the absolute value 



of the cross section, but also the size of the higher-order QCD corrections. Since the accompanying 
QCD radiation is softer at the Tevatron than at the LHC, the impact of the veto is less effective at 
the Tevatron than at the LHC. In the case of a strong transverse-momentum cut of p^ to = 15 GeV 
at the Tevatron, the NNLO contributions increase the NLO result by about 25%. At the LHC, the 
difference between the NNLO and NLO results is already reduced to less than 10% by considering 
a weaker cut of p™ to = 30 GeV. By further decreasing p™ to , perturbative contributions beyond 
the NNLO can become sizeable, since they are enhanced by powers of logarithmic terms of the 
type lnp™ to . 

In this paper we have limited our analysis to the SM Higgs boson. The NLO QCD corrections 
to Higgs boson production within the Minimal Supersymmetric extension of the Standard Model 
(MSSM) are known [jrj], |9j; for small values (tan/? ^5) of the MSSM parameter tan/3, they 
are comparable to those for the SM Higgs. This suggests that this similarity of the size of the 
perturbative QCD corrections remains true also at NNLO. 

In actual experiments at hadron colliders, the observed cross section for Higgs signal events 
can significantly depend on the details of the experimental analysis. This is particularly true when 
a jet veto procedure is applied, because of effects due to efficiency for jet reconstruction, jet energy 
calibration, presence of pile-up events and so forth. Quantitative estimates of these effects require 
event simulations H [j]] based on parton shower Monte Carlo generators f42"l , 43]. The parton 



shower produces multi-parton configurations that approximate the exact QCD matrix elements, 
but such approximation does not strictly correspond to the perturbative expansion of the cross 
section at LO, NLO, NNLO and so forth. Therefore, the NLO or NNLO fT-factors computed in 
this paper cannot naively be used to rescale the results of present Monte Carlo simulations at the 
the Tevatron and the LHC. As is customary practice in QCD analyses, more refined studies, which 
combine the perturbative QCD predictions with the Monte Carlo simulations at the detector level, 
are necessary to firmly quantify the expected number of veto selected Higgs events at Tevatron 
and the LHC. 
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